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Abstract —This research aims to study dynamic and static conditionsaspect deposition of particle stimulation 
applies Granular Dynamics model to describe contact forces simultaneously.The movement of particles forms a 
trajectory.The trajectory of particles resolved by Gear Predictor-Corrector method. Simulation involved a 
variation of thestandard deviation value of particle diversity in two dimensional simulation of polydisperse 
particles. The results of simulation areimages of deposition process for every time iteration until the pile of 
particles in static condition. Conclusion is the smaller standard deviation value, cause the greater angle of 
reposeof pile, and otherwise the greater standard deviation value, cause the smaller angle of repose of pile. So, 
the standard deviationsize diversity of particles is inversely proportional height and slope of forming pile 
structure. 
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I. INTRODUCTION 

Granular particles are related to daily life, such as sand, 
stones, grain of rice, beans, sugar, potatoes, and 
capsules.Granular particles behaviourare unique. They can 
be like a solid, liquid, and gas. It was involved with 
external condition [3].Particles dynamics are study of 
forces and motion that occur of particles. Forces of 
movement particles are normal force, tangential force, and 
gravity force. A lot of treatments from particle dynamics is 
deposition. Particles deposition are dropping grains into 
free surface, a pile is formed. The pile almost forms a 
triangular with irregular surface. 

Particles deposition has been studied force of collision 
model with stick and slip phenomenon. It result of multiple 
transition phenomenon at behaviour particle during 
deposition process. Multiple transition is from statics to 
dynamics and dynamics to statics. The simulation use 
monodisperse particle for simple model [6]. Not only 
monodisperse particle, but also we can simulated to find 
out the particle dynamics with different size, it call 
polydisperse particles. For polydisperse particles has 
discussed on [1] and [4].On reference [4], It has discussed 
polydisperse frictional particles packing. The conditional 
probability distributions of particle overlaps are 
determined by molecular dynamics simulations. It depends 
on polydispersity and friction coefficient. 

On the other hand, this deposition process involves 
many particles to be simulated. The particles which are 
simulated statistically are considered as data, which in this 
process can be observed the distribution of particles that 


occur. Discussion of particles within large numbers refers 
to analysis of particle size distribution. It is used to 
observe, control, and analize large quantities of industrial 
materials [3]. For another type of particles namely 
polydisperse hard sphere,On reference [1], It has discussed 
the comparison between particle size distribution with 
gauss distribution and result experiment of Transmision 
Electron Microcopy (TEM) method. They also did 
variation on standard deviations to determine of particles 
scattering. 

Based on previous researches, in this paper discuss 
polydisperse particles deposition with variation on 
standard deviations values. The purpose of writing is to 
find out the dynamic and static conditions result of 
dropping particle by providing variations in standard 
deviation. 

II. RESEARCH METHODS 

The first step of simulation was determining collision 
criteria. Collisions between particles occur when the 
distance between the two circles is less than the sum of the 
radius of the two particles. Formulation of collision criteria 
for polydisperse particles the number of the two radius is 
formulated with r x + r 2 because the size of radius particles 
is different. So, collision criteria for circular particles with 
heterogeneous diameter sizes formulated by (p t j = d t j — 
(fi + Tj). It given criteria cp t j < 0. d tj - is distance of center 
point of two particles, r t radius of particle i-th and rj radius 
of particle j-th. While the criteria of particle collision with 
media, for polydisperse particles cp zi formulated by cp zi = 
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z i —r i . It given criteria cp zi <O.Zi is the position of 
particle in vertical direction. 

The second step is simulation with Digital Visual 
Fortran software and interpretation of results. The number 
of particles simulated N = 900 particles by giving three 
different colors namely red, yellow, and green for each of 
the 300 particles. Particles dropped from height h = 3 m, 
normal and tangential velocity v n = v t = 5xl0 -3 ms -1 , it 
is assumed that all particles have the same mass m = 
0,05 kg althought the diameter is different, elasticity 
constants in the normal and tangential directions k n = 
k t = 10 s kgs~ 2 , stick velocity £ = 10 _2 ms _1 , time step 
At = 10 -5 s, coefficients of dynamics friction fi d = 0,3, 
and coefficients of statics friction fi s = 0,6. The value of 
normal damping constant y n has given by Equation 1 

Yn-crit — 2yjm. k n , (1) 

given parameter values as known previously y n - C rit ~ 
140. Because k n = k t , so 

Yt-crit — 2yjm. k t — Yn-crit' 

where Yn-crit critical normal damping constant and Yt-crit 
is critical tangential damping constant. Normal damping 
constant has a half value of y n _ crit , whereas tangential 
damping constant has equal with Yt-crit- Given y n = 0,5 
and y n - C rit equal with coefficient of restitution value (i)/) is 
0,3. Described i); in Equation 2 

^(*-2arctan@)) (2) 

1 

withoo = (4mk t - y 2 )2 [5]. 

The random number of diameter particles has obtained 
the largest of diameter particles i.e.0.121 m. Collision of 
two particles cause particles bounce to unexpected 
direction. So we can set the time collision of two particles 
to avoid it. For simple program, we use the largest 
particles for determine safest distance of particles. So, for 
the smaller particles are follow it. Formulated of mileage 
free fall particle is 

1 

s = -.g.t 2 + v 0 .t (3) 

where g is acceleration of gravity (9,81 m/s 2 ), t is time, 
and v 0 is initial velocityof particles. Ifs equal the largest of 
diameter particles and iff = n. AtwithAt = 10 -5 then 
based on Equation (3)obtain n « 15.706.We choose 
distances between first particle falling in order not to 
collide with the second particle must fall after 15.800-th 
iteration. The interpretation of the results dynamically is 
seen from a collection of images made by the movie with 
the movie maker software and the final image results in a 
pile in a static condition. 

Random particle diameters are generated using the 
Linear Congruent Generator (LGC) method, which is a 


randomly generated random number generation method 
using the gauss function. Normal distribution has two 
parameters namely mean p and standard deviation a. If we 
give the same mean but different standard deviations, 
curve are centered at exactly the same position on the 
horizontal axis, but the curve with the larger standard 
deviation is lower and spreads out farther [8]. Particle 
velocity is defined in the x-axis direction, whereas in the z- 
axis direction it is made zero so that all particles are given 
a velocity that has a direction. Although the diameter 
particle is random number, but it is assumed that all 
particle have the same mass. The purpose of that 
assumption is to focus of heterogeneity particles. 

Trajectory particles determinated by using the 
Predictor-Corrector method which is spesifically using the 
Gear Predictor-corrector method. As for the steps are 
prediction, evaluation, and correction stages [7]. Prediction 
stage is used to predict the initial position value of a 
particle, velocity, acceleration, third derivative, and fourth 
derivative of the particle position. Evaluation stage 
contains steps for calculating each force involved on all 
particles during the dynamic process, so that they can 
determine the possibility of particles colliding with each 
other. The Correction stage is the last step in Gear method. 
The acceleration of the correction is obtained from the 
force that has been obtained from the evaluation process 
divided by the mass of the particle, so that the values of 
position, velocity, acceleration, third derivative, and fourth 
derivative of the particle position are better. 

The force calculation starts with calculating the 
possibility of particles colliding with each other. Two 
colliding particles have their own criteria which have been 
explained previously. Other particles that are within a 
range of these particles are called neighbors, while 
particles that have the possibility to collide are called 
partners. The core program also defines particles that are 
more likely to collide, called neighbor list. It is efficient 
the running program by not needing to do calculations on 
n — 1 other particles. In addition to calculating the 
probability of collisions between particles, it is also 
examined the possibility of collisions of particles with the 
medium in which the particles are dropped. 

If particle qualify the collision criteria with a particle i 
then the particle will be save as a partner of particle i, then 
the normal force f n is calculated between the two particles 
by first calculating the relative velocity, relative 
acceleration, and the unit vector in the normal direction 
using Equation (4) 

fn ~ ~Q^n^n Yn^n)> (4) 

where 8 n is the relative displacement in the normal 
direction, v n is the relative velocity in the normal 
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direction, k n is the normal elastic constant, and y n is the 
normal damping constant. Then the particle displacement 
and relative velocity in the tangential direction are 
calculated to get the friction between the two particles 
using Equation (5) 


f = C^t^t + YtVt)> \v t \ < £ (static friction) 

1 { —v t g d \f n \, \v t \ > e (dynamic friction) ^ J 

where f t is tangential force, k t is the tangential spring 


constant, 8 t = v t (t) dx„ i.e. the total tangential 

displacement during time t — t 0 ,v t and v t is the relative 
tangential velocity and vector shape, p. d is the dynamic 
friction coefficient, y t is the tangential damping 
coefficient, and 8 is speed threshold [5]. 

Tangential and normal force are summed to get the 
resultant force received by the particle. This process takes 
place at any time dt until all particles dropped before 
particle i are examined against the collision criteria for 
particle i using Equation (6) 

/ = m5+^Cj, (6) 

where / is the contact force between particles, m is the 
mass of the particle, g is the gravitational force of the 
earth (9.81 ms~ 2 ), and c t is the contact force on the i-th 
particle [5]. If particle qualifycollision criteria with the 
medium, the normal and tangential forces of the particle 
are calculated like collision with other particles. These 
forces are added with the previous collision forces also. 
Another force involved on particles deposition is gravity. 
This force works on each particle andreduce force that 
leads to tangential direction. And so on until all particles 
have been dropped and all particles have remained or until 
they reach a predetermined step limit. 


III. RESULT AND DISCUSSION 
3.1 Dynamics and final structure of pile 

These observations describe the dynamic condition of 
particles deposition. It was used to produce movie, 
consisting of a seriesconcecutive figures of the pile 
configuration. Based on movie, the first particles dropped 
will hit the medium directly. Then followed by particles 
that fall afterwards, it will hit particles that have fallen first 
so that the particles accumulate with each other. 

The particles dropped earlier being covered by layers 
ofdrop later, and piling up. Furthermore, the particles 
dropped may not replace the earlier particle, but it go 
down through the edge of the pile directly and avalanches 
surface of the pile has been passed. After a long time, the 
pile is formed higher and parts of pile that avalanche are 
the edges. It is caused by the displacement of particles that 
fall later the pile is formed to obtain a shorter pile than the 
particles dropped earlier. So the force that occurs on 


particle cannot overthrow the pile has formed. Particle 
movement continues although all particles have been 
dropped. It caused by simultaneous movement of dropping 
particles, so that some of the particles are still in the 
neighbor list internal landslide until each particles at stable 
position.Result of simulation is stable pile of polydisperse 
particles deposition given by Fig.l 



Fig. 1 Stable pile of polydisperse particles deposition 


Based on Figure 1, small particles are more dominant 
in the middle of pile. This results also caused distribution 
of large-sized particles dominate the outside of pile, 
whereas the middle of pile is dominated by small-sized 
particles that push large particles into the outer side. The 
formed of pile is tends to be sloping, due to smaller sized 
particles entering between slits the larger sized particles. 

The pile can be estimated as isosceles triangle. It has an 
almost symmetrical triangular shape with an irregular 
surface. So, angle of repose the pile ^calculated by 
Equation (7) 



with h the height of triangle and b is the base of triangle. 
The height and base of triangle given by Equations below 
h = 3 y c 



y c = an d x c = x t . 

v z and y t are the coordinates of particles centre, x c and y ( are 
the centre of mass of the pile and the triangle [6]. So, the 
stable pile of polydisperse particles depositionas Figure 1 
was obtained angle of repose 27,45°. 

3.2Variation of standard deviation values 

Variation of standard deviation is given to find out the 
dynamic and static conditions result of dropping particles. 
It presented in Table 1. 
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Table 1. Comparison of the results simulation with 
variation standard deviation values 


Standar 
deviasi a 


Angle of 
Repose 

(°) 


Image 


1(T 3 



1(T 2 



10" 1 



Based on the Table 1, the smaller standard deviation 
value, cause the greater angle of repose of pile, and 
otherwise the greater standard deviation value, cause the 
smaller angle of repose of pile. So it can be concluded that 
the standard deviation values of particles diversity is 
inversely proportional to the angle of repose occurring on 
the pile. 

Associate results between simulation of deposition 
polydisperse particle and statistical approach with 
variation of standard deviation values, they have the same 
behavior. As a proof this statement, it given explanation 
with Figure3, difference between the stable pile of 
monodisperse and polydisperse particles deposition. 



(a) (b) 

Fig. 3 (a) Stable pile of monodisperse particles deposition, 
(b) Stable pile of polydisperse particles deposition 

The form of pile resultsdeposition monodisperse particles 
is higher than result of polydisperse particles. The angleof 
repose monodisperse particles is 52,85° and for 
polydisperse particle is 27,45°. 

IV. CONCLUSION 

The results of deposition polydisperse particle 
simulation are form of pile lower than 


monodisperseparticle deposition. The variation of standard 
deviation values in the program gives the fact that the 
physical concept of deposition of heterogeneous particles 
forms a pile that mathematically resembles the concept of 
a normal distribution curve. The standard deviation of the 
particle size diversity is inversely proportional to the 
height and slope of the formed pile structure. 
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